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Abstract 

Aims. We seek to understand the impact on shape estimators obtained from circular and elliptical shapelet models under two realistic conditions: 
(a) only a limited number of shapelet modes is available for the model, and (b) the intrinsic galactic shapes are not restricted to shapelet models. 
Methods. We create a set of simplistic simulations, in which the galactic shapes follow a Sersic profile. By varying the Sersic index and 
applied shear, we quantify the amount of bias on shear estimates which arises from insufficient modeling. Additional complications due to PSF 
convolution, pixelation and pixel noise are also discussed. 

Results. Steep and highly elliptical galaxy shapes cannot be accurately modeled within the circular shapelet basis system and are biased towards 
shallower and less elongated shapes. This problem can be cured partially by allowing elliptical basis functions, but for steep profiles elliptical 
shapelet models still depend critically on accurate eUipticity priors. As a result, shear estimates are typically biased low. Independently of the 
particular form of the estimator, the bias depends on the true intrinsic galaxy morphology, but also on the size and shape of the PSF. 
Conclusions. As long as the issues discussed here are not solved, the shapelet method cannot provide weak-lensing measurements with an 
accuracy demanded by upcoming missions and surveys, unless one can provide an accurate and reliable calibration, specific for the dataset 
under investigation. 

Key words, gravitational lensing - techniques: image processing 
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1. Introduction 

Shapelets have been proposed to describe the complex mor- 
phology of galaxies and to compress this information in a small 
number of expansion coefl&cients (Refregier 2003). Indeed, 
Massey et al. (2004) find that the shapelet decompositions 
of deep HST observations preserve the distribution of several 
morphological estimators so well that one can create simu- 
lations of a mock sky which are essentially indistinguishable 
from actual observations. 

For weak-lensing measurements it is furthermore important 
that shapelets can also faithfully describe the PSF shape (e.g. 
Bernstein & Jarvis 2002; Refregier 2003; Massey & Refregier 
2005; Kuijken 2006). Additionally, deconvolution from the 
PSF is an analytic operation in shapelet space (Refregier & 
Bacon 2003; Melchior et al. 2009). This means, given shapelet 
models of the observed galaxy and the PSF, respectively, only 
linear algebra is required to obtain a shapelet model of the de- 
convolved galaxy shape. Finally, several estimators for gravita- 
tional shear and flexion can be directly formed from this model 
(Bernstein & Jarvis 2002; Massey et al. 2007). In sunmiary, 
once a shapelet model of each galaxy is obtained, all necessary 
steps for weak-lensing measurements can be performed very 
eflficiently in shapelet space. 



There is, however, an implicit assumption in the setup of 
such a shapelet-based weak-lensing pipeline: namely that the 
occurring shapes can be faithfully described by shapelet mod- 
els. But this assumption is strictly valid only if the basis system 
is complete, which requires an infinite number of modes. In 
particular in weak lensing, the amount of modes which can be 
measured from the images is massively limited. For a limited 
number of modes, shapes which deviate from the characteris- 
tics of the employed basis functions can only have approximate 
representations. 

This limitation applies to any complete basis function set. 
In this work, we want to understand the consequences partic- 
ularly for weak-lensing measurements within the shapelet ba- 
sis. We start by describing the simulated galaxy images in sect. 
2. As there are two general shapelet variants, namely circu- 
lar and - as a generalization - elliptical ones, we discuss both 
of them in sects. 3 and 4. We show the influence of PSF con- 
volution, pixelation and pixel noise in sect. 5, and discuss the 
consequences of our findings in sect. 6. 



2. Test images 

For this work we seek to answer the following question: How 
well can a shapelet model recover the shear from a source 
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whose shape is not entirely describable by the model. To pin 
down the effects of an incomplete shape description and to 
isolate it from other systematics, we generate very simple test 
cases in which the galaxy shapes are - initially - not affected 
by PSF convolution and pixel noise. 

The radial light profile of galaxies is often described by the 
Sersic profile (Sersic 1963) 

p^(r) cc cxp{-K\{^f"' -l]], (1) 

where Rg is the radius containing half of the flux^ and Us is 
the so-called Sersic index. This type of profile is identical 
to a Gaussian for rig = 0.5 and steeper for ris > 0.5. From 
several investigations it is known that for the vast majority of 
observable galaxies Us ranges between 0.5 and 4 (e.g. Sargent 
et al. 2007)^. Thus, we describe the intrinsic shapes G with a 
flux-normalized Sersic profile with Re e {5, 10,20} pixels and 
Hs e {0.5, 1, 2, 3, 4}. For ensuring finite support, G is truncated 
Sit 5 Re. 

G is sheared in real-space by transforming the coordinates 
by means of the linearized lens equation (e.g. Bartelmann & 
Schneider 2001) 

"('-;r .7;,)-' 

such that G'(x') is obtained by evaluating G at the position x. 
The values of yi range between and 0.5; 72 is set to zero. It 
is important to notice that G is circular, while observed galax- 
ies show a wide distribution of intrinsic ellipticities (Bernstein 
& Jarvis 2002). Hence, G has to acquire its intrinsic ellipticity 
entirely from the applied shear. To obtain roughly realistic re- 
sults, the applied shear is varied up to I7I = 0.5 although such 
values cannot be generated by the cosmic large-scale structure 
and are even untypical for all but the innermost parts of galaxy 
clusters. As an advantage of this procedure, G has elliptical 
isophotes, for which the axis ratio and orientation are consis- 
tent at all radii, and therefore all ellipticity measures formed 
from these images should agree. 

G' is sampled at the final resolution of 20 Re x 20 Re pixels. 
While Re = 5 is already rather large for typical weak-lensing 
galaxies, we chose to simulate also even larger ones so as to 
mimic higher resolution images from which we can assess the 
impact of pixelation on the shear estimates. 

^ This is ensured by demanding of b^^ to satisfy the relation 

r(2n,) = 2y(2n,,bn,) 

between the complete and the incomplete Gamma function (Graham 
& Driver 2005). The approximate solution for the equation above, 

bn, ^ 1.9992^,-0.3271 

(Capaccioli 1989) is used throughout the paper. 

^ The exponential profile which describes the brightness distribu- 
tion of spiral galaxies and the de Vaucouleurs profile (de Vaucouleurs 
1948) of elliptical galaxies are special cases of the Sersic profile with 
ris = I and Us = 4, respectively. 



Because there is no pixel noise in these test images and 
the resolution is very high, the ellipticity e measured from 
quadrupole moments of the pixelated image is always compat- 
ible with the shear 7. Additionally, the centroid position can be 
computed with essentially arbitrary precision from the image. 

3. Circular shapelets 

According to Refregier (2003), 

Bn(x;/3) = /3-'(f>n,(fi-'xi) (f>n,(fi-'x2) (3) 

defines the two-dimensional shapelet basis function of order 
n = (^1,^2) and scale size j3 which is related to the one- 
dimensional Gauss-Hermite polynomial 

0„(x)^[2"7r^^!]-^if„(x)e-T, (4) 

with Hfi being the Hermite polynomial of order n. 

From the shapelet coefficients of some two-dimensional 
function G 

d^xGix)Br,{x;fi), (5) 

00 

one can reconstruct a shapelet model 

ni+n2-n,nax 

G(x)= £ gnB„ix;/3). (6) 

ni,n2=0 

The model is found by minimizing 

X^c^Yj^G(x)-G(x)f (7) 

X 

with respect to the parameters gn. JS and rimax could be fixed 
at values chosen to be suitable to describe the galaxy ensemble 
well on average (Kuijken 2006) or be determined from the min- 
imization of as well (Massey & Refregier 2005; Melchior 
et al. 2007). 

One foreseeable problem of the shapelet decomposition 
stems from the Gaussian weighting function in Eq. (4). We dis- 
cussed already that galaxies typically have steeper profiles than 
a Gaussian, which means that an optimized shapelet model re- 
quires higher orders. However, due to the polynomial in Eq. (4), 
the largest oscillation amplitudes of high-order modes are lo- 
cated at rather large distances from the centroid. Models which 
include higher orders thus allow a better description of the outer 
parts of a galaxy, while they still fail to reproduce correctly 
the central region in the case of steep profiles. Additionally, in 
case of noisy image data, the number of modes must be limited 
to avoid overfitting spurious nearby noise fluctuations. Hence, 
galactic shapes with steeper profiles than a Gaussian are ex- 
pected to be described by shapelet models with systematically 
shallower profiles. 

Of similar concern is the circularity of the shapelet basis 
system. As the scale size for both dimensions in Eq. (3) is the 
same, the zeroth-order is round. If the shape to be described is 
stretched in a particular direction - as a result of its intrinsic 
shape or due to gravitational lensing - this elongation has to be 
carried by higher shapelet orders. Again, for a limited number 
of basis modes we must expect an insufficient representation 
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of the true shape by the shapelet model. In particular, we have 
to consider an underestimation of the source elongation much 
more likely than an overestimation, as the basis system prefer- 
entially remains circular. 

These shortcomings are likely but do not have to affect the 
estimation of a weak shear. Whether a model with limited fi- 
delity leads to a biased shear estimate depends an the construc- 
tion of the estimator. 

The most straightforward estimator is obtained from the 
complex ellipticity (e.g. Bartelmann & Schneider 2001) 

Qn-Q22^2iQn 

where Qij are the quadrupole moments measured from the 
shapelet model & . From this we obtain a direct shear estima- 
tor 7^^^ = 6. By choosing this particular definition, we benefit 
from an ellipticity estimator with a perfect response to external 
shear, i.e. the shear responsivity R = de/dy (e.g. Bernstein & 
Jarvis 2002) equals unity and we do not have to correct for it. 

It is important to note that the Qij integrate over the en- 
tire shape of G\ consequently Qij are linear combinations of 
all available shapelet coefficients of & (Berge 2005), therefore 
this estimator critically relies on a decent shapelet model. 

3.1. Polar shapelets 

For the problem of describing galactic shapes, it is often more 
convenient to use a polar coordinate frame (r, (p) instead of the 
Cartesian frame (xi, X2). Following again Refregier (2003), we 
define the polar shapelet basis functions as 

Bn,,nM^ (f) = p-'(l)nr,miP~'r, if), (9) 

where ni and rir are the left-handed and right-handed modes, 
respectively, and 

(l>n.,nXr,<p) ^ Vnnr\ni\r^Hn,,nXr)e-^^e'^-^-''^^^. (10) 

Bernstein & Jarvis (2002) showed that, for ni < rir, one can 
relate 

Hn,nXr) ^ (-l)"'n/!r"-"'L«;-Xr') (11) 
to the associated Laguerre polynomial 

L^W^^^Cr^-^-O. (12) 

^ p\ drP 

Note the similarities between the Cartesian basis Eqs. (3) & 
(4) and the polar basis Eqs. (9) & (10): Both share a Gaussian 
weighting function and are intrinsically circular. 

As Bn^^m is a complex function if riy nu it is computa- 
tionally more efficient to decompose the galactic shape into 
Cartesian shapelets according to Eqs. (5) - (7) and then to per- 
form a coordinate transformation T^^^ from Cartesian to polar 
shapelet space (cf. Eq. (69) in Refregier 2003), 

8nr,ni ~ T ^ gn\,n2'> (1-^) 

We will do that for the tests in the following sections. 

Performing this transformation enables us to form a con- 
ceptionally diff'erent family of shear estimators. By defining 



n = rir + ni and m = riy - we can see from Eq. (10), that 
Bn^fn behaves like a field with spin m. Since the shear behaves 
like a spin-2 field, its action on a circular source can be de- 
scribed by the m = 2 modes of a polar shapelet model. Thus, 
the most basic one of these estimators, 

fp 

= 1 ^ (14) 

V^(^ + 2) <g^_2,o - C2,0> 

uses only the m = 2 mode of any even polar order n of & 
(Massey et al. 2007)^. This estimator must be normalized by 
radial modes g^^ obtained from unlensed sources. For the tests 
discussed here, we provide the unlensed image G such that 
those coefficients can be measured from images having the cor- 
rect unlensed shape. 

This approach is drastically different from measuring the 
quadrupole moment - in the sense that it depends only on a 
single lensed and two radial unlensed coefficients - and thus 
expected to behave in a different manner, even though the de- 
composition is done with Cartesian shapelets in both cases. 

3.2. Results 

G' is decomposed into Cartesian shapelets of maximum order 
J^max ^ {8, 12}, which is typical given the significance of weak- 
lensing images (cf. Kuijken 2006). 

At first, we investigate the modeling fidelity visually. In 
Fig. 1, we give four examples of Sersic-type galaxy shapes and 
their shapelet models. It is evident from the top row that for 
modest ellipticities, an elliptical Gaussian can be very well rep- 
resented by its shapelet model. But if either the ellipticity be- 
comes stronger or the intrinsic galactic profile becomes steeper, 
the shapelet decomposition performs poorer. For the Gaussian 
case shown in the third row, the overall shape is evidently more 
compact and rather boxy than elliptical, and aff'ected by oscil- 
latory artefacts. The images with Hs = 2 (second and forth row) 
show prominent ring-shaped artefacts and are concentrated at 
the core region of G. It is striking that the drastic increase in 
ellipticity from €\ = 0.15 to 6i = 0.5 causes almost no apparent 
change in the respective shapelet models. 

The same trend can be seen quantitatively in Fig. 2, where 
we show the logarithm of the goodness-of-fit;^^ of the shapelet 
models. From this we infer that the modeling errors become 
stronger with increasing ris and are almost independent of y 
for Hs > 1. In other words, while the shapelet decomposition 
is occupied with modeling a steep profile, it misses most of 
the ellipticity information. It is worth noting that an increase 
of Umax from 8 to 12 does not lead to substantially lower 
although the number of available modes is raised from 45 to 
91. This behavior can be explained by the shape of the higher- 
order shapelet functions: As they preferentially fit features in 
the outer regions, they improve the model only in the low-flux 
regions. 

In the top panel of Fig. 3, we show the shear estimator 7^^^ 
as a function of the applied shear. We see that the bias is es- 
sentially a linear function of yi with a negative slope which 

^ We restrict ourselves to n = 2 for the following tests, so that we 
only discuss y^^^\ 
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n, = 0.5, 61 = 0.15 



n, = 2, ei= 0.15 



: 0.5, 61 = 0.5 




Us = 2, €i = 0.5 






Figure 1. (left) Examples of Sersic-type galaxy images with 
Hs = 0.5 or 2 and an intrinsic ellipticity (induced by shearing 
the circular profile Eq. (1)) of 6i = 0.15 or 0.5. (right) Best-fit 
circular shapelet models & with rimax = 12. The color stretch 
is logarithmic. 



increases with n^. It is important to note that for Us = 0.5 
the estimator is unbiased as long as y remains moderate. But, 
for Us > 3 the same estimator is essentially shear-insensitive. 
Increasing the maximum order rimax from 8 to 12 improves the 
estimator, because the shape at large distances from the center 
is captured better by the model, and by construction f^^^ makes 
use of all available orders. But the steeper the profiles becomes 
the less can high orders contribute to the shear estimation, be- 
cause the quadrupole moments become dominated by the inner 
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Figure 2. Decomposition goodness-of-fit as function of the 
applied shear for Sersic-type galaxy images with 0.5 < < 4. 
The shapelet models use rimax = 8 (dashed lines) or 12 (solid 
lines). As the images are noise-free, the units of;^^ are arbitrary. 

region which is governed by a single central pixel with square 
shape and therefore vanishing ellipticity^. From our prior dis- 
cussion and Figs. 1 & 2, we anticipated such a behavior and 
it clearly confirms our theoretical expectations regarding steep 
galactic profiles. 

In the bottom panel of Fig. 3, the response of the shear es- 
timator 7^^^^ on the same set of galaxies is shown. We can see 
that the overall bias is mitigated by roughly a factor 4, rather 
independent of Umax- This has to be expected as the shapelet 
basis is orthogonal and thus higher-orders do not change the 
value of ^2^2' which carries the shear signal of f^'^^\ The diff'er- 
ences which occur when changing n^ax are related to a diff'erent 
preferred scale size JS in the optimization. 

As before, galaxies with ris = 0.5 can be measured with 
high fidelity. For steeper profiles the estimator has a - some- 
what surprising - positive bias, while the shapelet model itself 
underestimates the ellipticity. We do not fully understand why 
this estimator overestimates the applied shear, but we can iden- 
tify two possible reasons: 

- The estimator has been derived from the action of a in- 
finitesimal shear on a brightness distribution which is per- 
fectly described by a shapelet model (Massey et al. 2007). 
In the tests performed here, we intentionally violate these 
unrealistic assumptions. 

- Looking at the definition in Eq. (14) and the shapelet mod- 
els in Fig. 1, we see ring-shaped artefacts for steep pro- 
files, corresponding to radial shapelet modes. Exactly these 
modes are required for normalizing the estimator. As we 
know from Fig. 2, the goodness-of-fit - and hence the abun- 
dance of artefacts - is highly correlated with rig. Thus, the 
denominator of Eq. (14) is probably also plagued by the 
poor reconstruction quality of steep profiles. 

We can confirm the last argument by looking at the results of 
the next higher-order estimator 7^^^\ and find it to be strongly 
biased and highly unstable under variation of rig. This trend 



We will discuss the effects of pixelation in sect. 5.2. 
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^1 = n 

Figure 3. Shear estimates 7^^^ (top panel), 7^^^^ (bottom panel) 
for Sersic-type galaxy images as a function of the applied shear. 
The circular shapelet models use Umax = 8 (dashed lines) or 12 
(solid lines). 



continues for even higher polar order n and renders the family 
of estimators described by Eq. (14) unpredictable and thus un- 
usable for n > 2. Nevertheless, 7^^^^ is significantly less biased 
than7^^\ 



4. Elliptical shapelets 

In order to assess the performance of all versions of shapelet 
image analysis, we also considered an elliptical implementa- 
tion of the Gauss-Laguerre polynomials. To this purpose, we 
used a novel code that we have developed recently (Lombardi 
et al., in prep.). This implementation has many diff'erent op- 
erating modes regarding a number of key parameters, such as 
the truncation mode of the shapelet decomposition, the fitting 
of the object shape, and the masking method used. However, 
for sake of simplicity, here we mostly discuss the results ob- 
tained by running the pipeline in a mode that follows closely 
the prescriptions of Bernstein & Jarvis (2002) and Nakajima & 
Bernstein (2007). We refer to these two papers for all details 
that are not explicitly mentioned here. 

Elliptical shapelets are defined essentially like polar 
shapelets in a suitably sheared reference system. For example, 
in the simple case where the elliptical base is oriented along the 
Cartesian axes, one can define a new coordinate system (x^^x!^) 



as 

x[ = xi/a , X2 = X2/b , (15) 

where a and b are two scales (the two semi-axes of the base 
system). This gives rise to a radial coordinate for which we 
evaluate the shapelet function Bn^^m of Eq. (9) on this new sys- 
tem (by keeping JS = I, since JS is already encoded in the two 
scales a and b). Generalization for arbitrarily oriented galaxies 
is straightforward. 

All further steps in the decomposition are essentially iden- 
tical to the circular shapelet one: for example, the shapelet co- 
eflftcients are still evaluated using a technique identical to 
Eq. (7). 

The elliptical shapelet decomposition just described, and 
the consequent shape estimation, is done using a sequence of 
individual steps in our pipeline. 

1. First, we run SExtractor (Bertin & Arnouts 1996) on the in- 
put image (and associated weight map, if available) using a 
rather standard configuration. For each object, among other 
things, we measure its centroid, its shape, and its position 
angle. 

2. The resulting SExtractor parameters of each objects are 
used as initial guesses for the shapelet decomposition. We 
first extract the images of each object, and mask them us- 
ing an elliptical mask that follows the original SExtractor 
shape estimate (in particular, the mask ellipse has axes that 
are four times as big as the SExtractor ones). We then use 
the same SExtractor shape estimate as an initial guess for 
the shapelet elliptical basis. 

3. We perform a forward fitting of all pixels contained within 
the mask, and obtain this way a set of shapelet parameters 
g'n,m and associated covariance matrix. To be able to com- 
pare with the circular shapelets, we used a triangular trun- 
cation scheme by requiring that n = rir -\- m < rimax, where 
^max - 8 or 12 is the maximum order of the decomposition. 

4. The next step is the estimation of the ellipticity. For this 
purpose, we can either use the quadrupole moments as in 
Eq. (8) or rely on an iterative approach that we call "focus- 
ing": we apply a sequence of translations, scaling opera- 
tions, and shear operations to the object until it appears per- 
fectly centered, with maximum signal to noise, and round. 
As shown by Bernstein & Jarvis (2002), these conditions 
can be fulfilled by requiring that g^^ = g^^ = g^^ = 0. 
Since we are using an elliptical base for our shapelets, the 
condition g'^^ = simply indicates that the object is round 
in the elliptical base, i.e. has an ellipticity and a position 
angle that are identical to the ones of the base. 

5. Finally, we compare the focusing parameters with the orig- 
inal SExtractor one, and we repeat the whole process with 
the newly determined parameters in cases both parameter 
sets diff'er significantly. This step is typically executed for a 
relatively small fraction of galaxies only. 

A key point of the whole process is the presumed indepen- 
dence of the final result on the input SExtractor parameters. Of 
course, since the ellipticity estimation involves two non-linear 
iterative approaches (one for the design of the shapelet base 
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0.0 0.1 0.2 0.3 0.4 0.5 

^1 = n 

Figure 4. Shear estimate 7^^^^ for Sersic-type galaxy images as 
a function of the applied shear. The elliptical shapelet models 
use Umax = 8 (top panel) or 12 (bottom panel). The dotted line 
in both panel represents the bias of -20%, which we artificially 
applied to the initial guess of the ellipticity. 

and image mask, and a second for focusing), we cannot ex- 
pect a perfect independence. In other words, if we start from 
completely wrong source centroids and shapes, it will be very 
unlikely that we will be able to perform decent decomposi- 
tions and that the whole pipeline will converge towards sen- 
sible results. On the other hand, if the method relies too much 
on "good" initial parameters, and does not tolerate relatively 
small errors on the SExtractor catalog, we will likely face se- 
vere problems, because although SExtractor has a very good 
performance in common cases, it clearly has not been designed 
with weak lensing studies in mind, and does not reach the ac- 
curacy needed in this field. 

4.1. Results 

In order to test the ability of our elliptical shapelet pipeline to 
deal with inaccurate input catalogs, we biased the input elliptic- 
ity by 20% toward circular objects, and measured the residual 
bias left in the recovered ellipticity y^^^\ 

As shown by Fig. 4, this test produced conceptually simi- 
lar results as with circular shapelets: for small Sersic indices, 
we could recover the true ellipticity without any significant 
bias, while the estimates degrade significantly as we approach 
Hs = 4. The situation is improved when rimax is raised, because 
the transformations done during the focusing step take all avail- 



able orders into account, in contrast to the simpler description 
underlying the construction of y^"^^ in Eq. (14). With rimax = 12 
(bottom panel of Fig. 4), galaxies with ris = 0.5 and 1 have 
shear estimates without bias. When raising Us beyond that, the 
bias is at first positive before it becomes negative. Investigating 
this feature more closely, we find this estimator to have an os- 
cillatory tendency once the bias sets in. 

In summary, it appears that the pipeline is unable to fully 
correct inaccuracies in the SExtractor catalog for galaxies 
showing a profile close to the De Vaucouleurs one. 

5. Observational systematics 

In more realistic simulations or observational data, the galactic 
shapes are recorded after convolution with the PSF, pixelation 
by the CCD, and degraded by pixel noise. We now discuss the 
impact of these eff'ects on shear estimation with shapelets. 

5.1. PSF convolution 

Clearly, a convolution creates shallower profiles which can 
be better described by shapelet models. Therefore, the typi- 
cal goodness-of-fit values, in particular for steeper profiles, are 
considerable lower than in the unconvolved case. If the PSF 
shape is perfectly described by its shapelet model, one can 
exactly undo a convolution in shapelet space. In such a case, 
the shape obtained by deconvolving a PSF-convolved galaxy 
model must approximate the true, unconvolved shape G' better 
than its direct model 

To verify this new hypothesis, we convolved Sersic-type 
galaxies in pixel space with PSF shapes P obtained from 
shapelet models^, 

C = PirG\ (16) 

For circular shapelets, C is modeled with shapelets and explic- 
itly deconvolved from P in shapelet space, while for the ellip- 
tical shapelets we obtain the unconvolved shape by convolving 
the model with the PSF and fitting the outcome to the image 
data. 

In Fig. 5 we compare the bias of the shear estimates from 
unconvolved galaxies images and from the same set of galax- 
ies after convolution with a Gaussian PSF with a FWHM of 5 
pixels - thus FWHM = Re. As expected, the overall shape of 
the convolved images - and thus also their unconvolved shapes 
- can be better modelled with shapelets. Hence, the shear esti- 
mator 7^^^ profits manifestly because it makes use of the entire 
shape information: its bias is lowered by a factor 3. We ob- 
served similar improvements also for y^^^^ and y^^^\ 

We varied the size of the PSF and found that in general 
the bias is lower for larger PSFs. Additionally, we found y^^^^ 
and 7^^^^ to be rather sensitive to changes of rig and FWHM/7?g. 
Finally, the bias of all estimators also depends in a non-trivial 
way on the shape of the PSF, because the fidelity of the model 

^ For our argument, we only exploit that convolution with the PSF 
renders observed profiles shallower. Hence, the last statement is prob- 
ably still true for imperfect shapelet models of the PSF - which are 
likely to introduce other systematics. 

^ In the terminology of this work, that means P - P. 
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^1 = n 



Figures. Similar to Fig. 3: Solid lines show bias on shear es- 
timate 7^^^ from circular shapelet models with rimax = 12 for 
unconvolved galaxy images, while dashed lines are obtained 
for the same set of galaxies after convolution with a Gaussian 
of 5 pixels FWHM. 



for the convolved galaxy image determines the accuracy of the 
estimators. 

5.2. Pixelation 

Images from CCDs are obtained by collecting the light within 
pixels of approximately square shape. For measuring shapes, 
pixelation has important consequences. If the size of the ob- 
ject is small compared to the pixel size, we cannot describe 
the true continuous shape of the object but rather its piecewise 
approximation with pixel- sized step functions. Modeling ap- 
proaches like the shapelets method can take pixelation into ac- 
count by integrating the model values within the pixels (Massey 
& Refregier 2005). In case of convolved images, the deconvo- 
lution procedures also treat pixelation consistently, if the PSF 
shape has been measured from images with the same pixelation 
(e.g. Bridle et al. 2008). 

For estimating the shear, an additional problem is of relevance. 
As the smallest piece of information within an image is given 
by a single pixel with square shape, we can only infer shear 
information from an object for which we can measure more 
than a single pixel. Particularly for galaxies with steep profiles, 
the largest fraction of the flux is registered in the pixel which is 
closest to the centroid. Then, also the shear information is dom- 
inated by this central pixel, which does not have any preferred 
direction, and hence biased low. 

For this work, it is important to verify that the biases re- 
lated to steep profiles are not entirely a pixelation problem, but 
stem from the shape mismatch as discussed before. Therefore, 
we also made sets of images with Re = 10, 20. 
Although there are some diff'erences between the three tested 
estimators, we found a common trend when increasing the size 
of the galaxies: The results for galaxies with ris < I are essen- 
tially unchanged. In particular, the bias does not vanish when 
the side length of a pixel is reduced to one quarter. This shows 
that there is a profile dependent bias, independent of pixelation. 



The estimates for steeper profiles benefit from smaller pixel 
sizes, but the results still remain more strongly biased than for 
shallower profiles. 

5.3. Pixel noise 

There is an additional eff'ect related to the discussion of pixela- 
tion: In the presence of pixel noise, a smaller number of signifi- 
cant pixels remains for each galaxy. In particular, steep profiles 
therefore tend to be reduced to some pixels or even only a sin- 
gle pixel close to center of the galaxy, which is then fitted by 
the model. Thus, we expect pixel noise to behave in a similar 
way as strong pixelation. 

We performed the same tests again with realistic noise 
added to G\ but apart from additional statistical uncertainty 
it did not lead to diff'erent results. In fact, we can confirm that 
steep profiles are aff'ected more strongly by pixel noise than 
shallower ones. 

There is another important point to note: Increased pixel 
noise would normally lead to a lower rimax^ as the shapelet mod- 
els are typically tuned such that they do not or only marginally 
pick up noise fluctuations. Reducing the maximum order typ- 
ically leads to more prominent modeling problems, as we al- 
ready discussed before and consequently to poorer results from 
most shear estimators. 

6. Conclusions 

- Shear estimates from circular shapelets are biased if the 
shape to be described has too steep a profile (steeper than a 
Gaussian) or too large an ellipticity. Profile mismatch is the 
more important source of bias. 

- For elliptical shapelets, profile mismatch still poses a con- 
siderable problem, because the shapelet models cannot 
fully correct biases of the ellipticity prior when the profile 
becomes steeper than exponential. 

- Diff'erent shear estimators can mitigate the bias, but never 
eliminate it completely because the shape mismatch gener- 
ally aff'ects all shapelet modes. 

- For flexion estimates we are even more concerned. While 
the overall octupole power might be lower than the 
quadrupole power, the bias introduced by the shape mis- 
match becomes more important for higher-order moments. 

- Convolution with a PSF renders all observable shapes shal- 
lower and allows the treatment of pixelation, and hence fa- 
cilitates a more accurate description by shapelet models. 

- The ellipticity estimates could be calibrated provided we 
have simulations detailed enough to properly describe the 
distribution of intrinsic galaxy morphologies and the efl'ect 
of PSF convolution on it. Furthermore, we would need to 
be able to infer the appropriate calibration parameters from 
observations. To achieve the latter is at least as hard as mea- 
suring the ellipticity. 

- Even if an efl'ective calibration for an entire galaxy en- 
semble could be achieved, it does not ensure unbiased re- 
sults when the ensemble is sliced (e.g. in redshift bins for 
weak-lensing tomography or spatially for cluster-mass re- 
constructions). 
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- Lensing simulations with shapelet source models (e.g. 
Massey et al. 2004; Meneghetti et al. 2008) systematically 
favor shapelet pipelines for weak-lensing measurements as 
the shapelet basis forms the optimal way of shape descrip- 
tion. 

- The mentioned shortcomings are mostly relevant for small 
and faint galaxies which are described by shapelet mod- 
els with low Uynax' For objects with a high signal-to- 
noise ratio, modeling with high Umax reduces these effects. 
Additionally, the artefacts which can occur on bright galax- 
ies have a similar impact on shapelet models for galaxies 
which are intrinsically similar. We conclude therefore that 
morphological investigations of galaxy populations are still 
feasible with shapelets. However, measures of e.g. their el- 
lipticity or profile slope distribution are likely biased. 

- PSF shapes typically follow the Moff'at profile (Moff'at 
1969) pmOc{\+ar^)-^- mid can thus be very well modeled 
already with a pure Gaussian - which is the limiting case 
for ^ oo (Trujillo et al. 2001). In particular for PSFs 
from ground-based observations with its at most mild ellip- 
ticity, shapelets can reach excellent modeling accuracy. 
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